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Abstract 

We theoretically investigate the effects of charge order and spin frustration on the spin ordering in TMTTF salts. 
Using first-principles band calculations, we find that a diagonal inter-chain transfer integral tq\, which causes spin 
frustration between the inter-chain dimers in the dimer-Mott insulating state, strongly depends on the choice of 
anion. Within the numerical Lanczos exact diagonalization method, we show that the ferroelectric charge order 
changes the role of tqi from the spin frustration to the enhancement of the two-dimensionality in spin sector. The 
results indicate that tq\ assists the cooperative behavior between charge order and antiferromagnetic state observed 
in TMTTFaSbFg. 
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1. Introduction 

Low-dimensional molecular conductors provide a 
fruitful stage to study strong electron correlation ef- 
fects leading to a wide variety of phase transitions [1]. 
In this context, among the most-studied families are 
the quasi-one dimensional (QID) molecular conductors 
TMTTF2X (TMTTF: tetramethyl-tetrathiofulvalene, 
X: monovalent anion) [1]. These salts form a QID tt- 
band at quarter-filling in terms of holes with intrinsic 
dimorization along the conduction axis. They exhibit 
various types of phase transitions such as ferroelectric 
charge ordering (FCO), spin-Peierls (SP), antiferro- 
magnetic (AF), and superconducting (SC) transitions 
by applying pressure or replacement of X [2-5] . Among 
them, (TMTTF)2SbF6 shows a pecuhar behavior un- 
der pressure; a cooperative reduction of FCO and AF 
phase transition temperatures by the application of 
pressure has been reported by NMR measurements [3] . 



This result naively does not coincide with the case for 
typical CO transitions, where CO suppresses the ten- 
dency toward magnetic ordering due to decrease of the 
effective spin exchange couplings [6,7]. Recently we 
have proposed a mechanism where the FCO transition 
increases the two-dimensionality in the magnetic sec- 
tor and leads to the stabilization of the AF state [8] . In 
this study, we have also found that the diagonal inter- 
chain hopping tq-i strongly depends on the choice of 
anion. Intriguingiy, tqi causes geometrical frustration 
between the chains of spins localized on every dimer, 
forming QID exchange couplings as seen from Fig. 1 
(b); we will see this effect more explicitly later. 

In this paper, we investigate the role of tq\ for 
the spin ordering in the absence or the presence of 
FCO, and how it affects the cooperative behavior be- 
tween FCO and two-dimensional (2D) AF states in 
TMTTF2SbF6. 
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Fig. 1. The schematic representation of the structures 
in the donor plane for (a) TMTTF compounds and (b) 
dimers composed of two TMTTF molecules. In dimer- 
ized pictures, the transfer integrals are effectively given by 
^ai — t'l^ — tb — tq2/2, and t'^-^^ — tqi/2, respectively. 

2 . For mulat ion 

We investigate a QID extended Hubbard model at 
1/4- filling in terms of holes with Coulomb interactions 
and inter-chain hoppings, whose Hamiltonian is given 
by 

■Hehm = - ^ tij{cl^Cj„ + H.c.) 

+U^^ni^nii + ^^Vijninj, (1) 

where tij is the transfer integral between the neigh- 
boring sites denoted by (ij), c]^ (dcr) is the creation 
(annihilation) operator of a hole on the ith site with 
spin a =t or J,, and rii — + with riicr = cl^Cia-. 
U and Vij are the on-site and the inter-site Coulomb 
interactions, respectively. 

The transfer integrals have been estimated [8] for 
two representative members, (TMTTF)2PF6 and 
(TMTTF)2SbF6 (hereafter, we call them as PFe salt 
and SbFg salt) by the first-principles band calcula- 
tions using the computational code QMAS (Quan- 
tum MAterials Simulator) [9] based on the projector 
augmented-wave method [10] with the generalized 
gradient approximation [11]. We determine the values 
by fitting to the band structures. 

The values of transfer integrals are given as 
{tai,ta2,tb, tqi,tg2} ={"155, -203, 26.2, -1.31, -3.29} 
for PFe salt and {-149, -207, 16.4, -16.4, -9.73} for 
SbFs salt, in the unit of meV (notations are shown in 
Fig. 1 (a) ). From these results, we can see that the di- 
agonal hoppings tqi,tq2 are negligible for PFg salt but 
appreciably large for SbFg salt, especially [tgi/tjl ~ 1. 

Considering these results, we hereafter set the trans- 
fer integrals as tai — —0.8, ta2 = —1, tt — 0.15, tq2 = 
and choose the diagonal inter-chain transfer integral 
tqi as a parameter. We set the on-site Coulomb inter- 
action as f/ = 4 and impose a constraint on Vij as 
Vai = Va2 ~ Vqi = Vq2 = V SLud Vj, = to realize the 
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Fig. 2. (a) Inter-chain real space spin-spin correlation functions 
between 1 and i sites (mirrti) (the site numbers are shown in 
Fig. 1 (a)) a.t V = for t,i = 0, - 0.15, and -0.3. Lines 
are connected between sites in a same chain, (b) A schematic 
figure for 2D AF state between dimers. (c) Intcr-dimcr spin 
structure factors S+(-7r, Qy) at V — for tgi — 0, — 0.15, and 
— 0.3. Lines arc guides for eyes. 

FCO pattern shown in Fig. 3 (b). 

We perform numerical exact diagonalization on a 
4x4 sites cluster under the periodic boundary con- 
dition and calculate the following inter-dimer/intra- 
dimer charge and spin structure factors given by 

C±iQ) = ^^J2i^My'''^'''~''''' (2) 

where A'^^ is the total number of dimers and denotes 
the center position of the ith dimer. Here, the inter- 
dimer/intra-dimer {-\-/ —) correlations are detected by 
the summation/difference in charge and spin densities 
within each dimer, — {n2i ± 7i2i-n)/2 and — 
{m2i ± m2i+i)/2 with mi — riif — nn, respectively, 
where even (odd) number is labeled as the site for the 
left (right) side in a dimer, as shown in Fig. 1 (a). 



3. Results 

First, we show the inter-chain real space spin-spin 
correlations {mirrii) at V = in Fig. 2(a). The sys- 
tem is in the dimer-Mott insulating regime [8] , where 
the intra-chain dimerization and U lead to a charge- 
uniform Mott insulator. For tqi = 0, the sign of {mirrii) 
changes each other between dimers; there exists AF 
correlation as shown in Fig. 2(b), although the inter- 
chain correlation is weak due to small ti,. With increas- 
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Fig. 3. (a) Intcr-chain real space spin-spin correlation functions 
between 1 and i sites {mirrii) (the site numbers are shown in 
Fig. 1 (a)) at y = 2 for t<,i = 0, - 0.15, and -0.3. Lines are 
connected between sites in a same chain, (b) A schematic figure 
for 2D AF state in FCO state. The size of a circle represents the 
hole density, (c) Intcr-dimer spin structure factors ^^(Tr, qy) at 
V — 2 for tqi — 0, — 0.15, and —0.3. Lines are guides for eyes. 

ing \tqi\, this AF correlation is diminished, as seen, 
for example, in the sign change of (mimg), for which 
the sites are connected by tqi. As a result, the AF 
correlation between inter-chain dimers is weakened. 
In Fig. 2(c), we show the inter-dimer spin structure 
factors S+{TT,qy). It has peaks at qy = ±7r, indicat- 
ing AF correlation between dimers. With increasing 
\tqi\, S+{tv, ±7r) are suppressed, while S+(-k, ±7r/2) and 
S'-i-(7r, 0) are enhanced; tqi operates as the inter-chain 
spin frustration parameter. 

Next, in Fig. 3(a), we show {mirrii) atV = 2 for the 
CO region. For tqi — 0, (m-imi) is large (note the dif- 
ference in the scale from Fig. 2(a)) only at odd num- 
bers iodd and its sign changes each other between the 
neighboring charge rich sites, indicating the enhanced 
development of 2D AF correlations in the presence of 
FCO, as shown in Fig. 3(b) [8]. When tqi is introduced, 
the absolute value of {mirrii^^^) becomes appreciably 
larger, in contrast to the case oiV = Q. Such a behav- 
ior is also manifested in S+{Tv,qy), as shown in Fig. 3 
(c), where S+{tv, ±7r) is seen to become large with in- 
creasing \tqi\, while S'+(7r, ±7r/2) and 5'+(7r,0) become 
small. These results indicate that tqi enhances the two- 
dimensionality in magnetic sector and further stabi- 
lizes the 2D AF state. 

The role of tqi on the 2D AF correlation depends 
on the degree of CO correlation. In Fig. 4, C_(0,0) 
and 5-i-(7r,7r) are plotted as a function of V for tqi — 
0, — 0.15, and —0.3. For each tqi, we see that C^{0, 0) 
and S'+(7r,7r) cooperatively develops in V > 1.5; the 
FCO correlation assists in stabilizing the AF state. The 
noticeable point is that tqi becomes relevant for the 
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Fig. 4. (a) Intra-dimer charge structure factors C_(0,0) and 
(b) inter-dimer spin structure factors S-|-(7r,7r) as functions of 
V for tqi — 0, — 0.15, and —0.3. Lines are guides for eyes. 

AF correlation in the FCO state, as seen from the en- 
hancement of S+ (tt, tt) by tqi only for V > 1.5, whereas 
the change is small for V < 1.5. On the other hand, tqi 
is irrelevant for the FCO correlation, i.e. C_(0, 0) does 
not depend on tqi . 



4. Discussion and Summary 

For insight into the effects of tqi on the magnetic 
properties, we consider the strong coupling limits and 
estimate the leading terms of the spin exchange cou- 
plings by perturbation calculations with respect to the 
transfer integrals. In the dimer-Mott insulating state 
for the case of large U and small V, the spin-exchange 
couplings between dimers along a-axis and 6-axis are 
given by Ja = —t'ii/Ud and Jb = —itl/Ud with the ef- 
fective on-site Coulomb interaction for the dimer units 
Ud [6]. Also, in the presence of tqi there is a diago- 
nal inter-chain spin-exchange coupling Jab = —tqi/Ud- 
This term causes the spin frustration between inter- 
chain dimers as seen from Fig. 1 (b) and is enlarged to- 
ward 2t(, ~ |iqi I; in fact, the qty dependence of 5+ (tt, Qy) 
becomes weak with increasing \tqi \ as shown in Fig. 2 
(c). 

On the other hand, in the basis of the FCO state 
for the case of large U and V, the coupling between 
the nearest-neighbor intra-chain charge rich sites is 
given by Ja ~ -4:tlitl2/{QUV^) from the forth-order 
perturbation [12]. For tqi = 0, the coupling between 
inter-chain charge rich sites is given by Jb ~ —At1/U 
from the second-order perturbation. Since Ja is sup- 
pressed with increasing V , Jb can become the same 
order compared to J a- This gives the increase of the 
two-dimensionality in the magnetic sector and induces 
2D AF state [8]. In the presence of tqi, there is another 
additional contribution for the coupling between inter- 
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chain charge rich sites from the third-order pertur- 
bation process (for example, 1 — >■ 8 ^ 5 — > 1 and 
1^5^8— >-lin Fig. 1 (a)), which is given by ~ 
—8tqitaitb/{SUV). Here, the minus sign comes from 
the minus term coupling with the transfer integrals in 
Eq. (1). Since tai is large compared to ti,, this process 
can become the same order compared to Jt in the in- 
termediate V region. As a result, this contribution be- 
comes relevant in the presence of FCO and the two- 
dimensionality in the spin sector increases. This agrees 
with the behavior of the V dependence of 5'+(7r,7r) 
shown in Fig. 4 (b). We can conclude that not only 
the enhancement of the two-dimensionality in spin sec- 
tor by FCO but also the largeness of tqi in SbFe salt 
compared to other salts promotes stabilization of the 
coexisting AF and FCO state observed in the experi- 
ments [3]. 

Finally, we comment on the role of tqi in the pres- 
ence of electron-phonon couplings, which is relevant 
for the SP state in both the dimer-Mott and the CO 
states [7,8,13,14]. Previous works show that, in general, 
the SP state is suppressed and the 2D AF state is sta- 
bilized with increasing the inter-chain spin exchange 
coupling Jt [15]. In PFe salt, it is experimentally in- 
dicated that, at finite-temperatures above the first or- 
der boundary between the two states, a quantum crit- 
ical behavior is seen, controlled by a hidden quantum 
critical point [16]. The competition between the SP 
and the 2D AF states occurs as well in the presence 
of FCO [3,8]. The role of tqi may enhance the strong 
quantum fiuctuations and bring about difficulties in 
determining the magnetic phases in the FCO state at 
low temperatures. One possibility is that the fluctua- 
tion leads to the destabilization of ordered phases in the 
critical region, and results in a spin-liquid state down 
to the ground state, in between the SP and the 2D AF 
quantum critical points. This tqi effect on the magnetic 
properties with electron-phonon couplings remains as 
a future problem. 

In summary, we have investigated the efi'ects of 
charge order and the diagonal inter-chain transfer 
integral on the spin ordering in TMTTF2X. The 
diagonal tqi works as spin frustration between the 
inter-chain dimers in the dimer-Mott insulating state 
and reduces the AF correlation. On the other hand, in 
the presence of FCO, this spin frustration is released 
and 2D AF correlation is developed; tqi enhances the 
two-dimensionality in the spin sector. This assistance 
of tqi for the cooperative behavior between the FCO 
and the AF states agrees with the experiments in 
TMTTFaSbFs. 
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